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Abstract — The availability of positional information is of great 
importance in many commercial, governmental, and military 
applications. Localization is commonly accomplished through the 
use of radio communication between mobile devices (agents) and 
fixed infrastructure (anchors). However, precise determination 
of agent positions is a challenging task, especially in harsh envi- 
ronments due to radio blockage or limited anchor deployment. 
In these situations, cooperation among agents can significantly 
improve localization accuracy and reduce localization outage 
probabilities. A general framework of analyzing the fundamental 
limits of wideband localization has been developed in Part I 
of the paper. Here, we build on this framework and establish 
the fundamental limits of wideband cooperative location-aware 
networks. Our analysis is based on the waveforms received at 
the nodes, in conjunction with Fisher information inequality. 
We provide a geometrical interpretation of equivalent Fisher 
information for cooperative networks. This approach allows us 
to succinctly derive fundamental performance limits and their 
scaling behaviors, and to treat anchors and agents in a unified 
way from the perspective of localization accuracy. Our results 
yield important insights into how and when cooperation is 
beneficial. 

Index Terms — Cooperative localization, Cramer-Rao bound 
(CRB), equivalent Fisher information (EFI), information inequal- 
ity, ranging information (RI), squared position error bound 
(SPEB). 



I. Introduction 

The availability of absolute or relative positional informa- 
tion is of great importance in many applications, such as 
localization services in cellular networks, search-and-rescue 
operations, asset tracking, blue force tracking, vehicle rout- 
ing, and intruder detection |[T|-|[8l. Location-aware networks 
generally consist of two kinds of nodes: anchors and agents 
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Fig. 1. Cooperative localization: the anchors (A, B, C, and D) communicate 
with the agents (1 and 2). agent 1 is not in the communication/ranging range 
of anchor C and D, while agent 2 is not in the communication/ranging range 
of anchor A and B. Neither agents can trilaterate its position based solely on 
the infonnation from its neighboring anchors. However, cooperation between 
agent 1 and 2 enables both agents to be localized. 



(see Fig. [T), where anchors have known positions while agents 
have unknown positions. Conventionally, each agent localizes 
itself based on range measurements from at least three dis- 
tinct anchors (in two-dimensional localization). Two common 
examples include the global positioning system (GPS) 191, 
ifTOl and beacon localization ITTI . lfT2l . In GPS, an agent 
can determine its location based on the signals received from 
a constellation of GPS satellites. However, GPS does not 
operate well in harsh environments, such as indoors or in 
urban canyons, since the signals cannot propagate through 
obstacles EJ-ID. Beacon localization, on the other hand, relies 
on terrestrial anchors, such as WiFi access points or GSM base 
stations ifTTl . lfT2l . However, in areas where network coverage 
is sparse, e.g., in emergency situations, localization errors can 
be unacceptably large. 

Conventionally, high-accuracy localization can only be 
achieved using high-power anchors or a high-density anchor 
deployment, both of which are cost-prohibitive and impractical 
in realistic settings. Hence, there is a need for localization 
systems that can achieve high accuracy in harsh environments 
with limited infrastructure requirements ||6l-|[8l. A practical 
way to address this need is through a combination of wideband 
transmission and cooperative localization. The fine delay 
resolution and robustness of wide bandwidth or ultra-wide 
bandwidth (UWB) transmission enable accurate and reliable 
range (distance) measurements in harsh environments ifTSl - 
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IfTSl Fl Hence, these transmission techniques are particularly 
well-suited for localization. Cooperative localization is an 
emerging paradigm that circumvents the needs for high-power, 
high-density anchor deployment, and offers additional local- 
ization accuracy by enabling the agents to help each other 
in estimating their positions H, IS, EH-lISl. In Fig. [T] for 
example, since each agent is in the communication/ranging 
range of only two anchors, neither agents can trilaterate its 
position based solely on the information from its neighboring 
anchors. However, cooperation enables both agents to be 
localized. 

Understanding the fundamental limits of localization is cru- 
cial not only for providing a performance benchmark but also 
for guiding the deployment and operation of location-aware 
networks. Localization accuracy is fundamentally limited due 
to random phenomena such as noise, fading, shadowing, and 
multipath propagation. The impact of these phenomena has 
been investigated for non-cooperative localization l?), ID, 
ll29l - ll3n . However, little is known regarding the bounds for 
cooperative localization. In particular, bounds on the coopera- 
tive localization performance were previously derived in IZTl . 
1281 using only specific ranging models. In other words, these 
works start from signal metrics, extracted from the received 
waveforms0 Such a process may discard information relevant 
for localization. Furthermore, the statistical models for those 
signal metrics depend heavily on the measurement processes. 
For instance, the ranging error of the time-of-arrival (TOA) 
metric is commonly modeled as additive Gaussian ETl . 1281 . 
13T1 . However, other studies (both theoretical IBl. l38l. l39l 
and experimental fS), ITSl ) indicate that the ranging error is 
not Gaussian. Hence, when deriving the fundamental limits of 
localization accuracy, it is important to start from the received 
waveforms rather than from signal metrics extracted from 
those waveforms. 

In Part I |[29l . we have developed a general framework to 
characterize the localization accuracy of a given agent. In this 
paper, we build on the framework and determine fundamental 
properties of cooperative location-aware networks employing 
wideband transmission. The main contributions of this paper 
are as follows: 

• We derive the fundamental limits of localization accuracy 
for wideband wireless cooperative networks in terms of 
a performance measure called the squared position error 
bound (SPEB). 

• We employ the notion of equivalent Fisher information 
(EFI) to derive the network localization information, and 
show that this information can be decomposed into basic 
building blocks associated with every pair of the nodes, 
called the ranging information (RI). 

• We quantify the contribution of the a priori knowledge of 
the channel parameters and the agents' positions to the 
network localization information, and show that agents 
and anchors can be treated in a unified way: anchors are 
special agents with infinite a priori position knowledge. 

'other aspects of UWB technology can be found in 1191 1251 . 

^Commonly used signal metrics include time- of-an i val (TOA) jvj, (§], 
flSl . fm . l32l . time-difference-of-arrival (TDOA) (33], |34|, angle-of-arrival 
(AOA) (3, EH, and received signal strength (RSS) Q, (36), E]. 



• We put forth a geometric interpretation of the EFI matrix 
(EFIM) using eigen-decomposition, providing insights 
into the network localization problem. 
> We derive scaling laws for the SPEB for both dense 
and extended location-aware networks, characterizing the 
behavior of cooperative location-aware networks in an 
asymptotic regime. 
The proposed framework generalizes the existing work on non- 
cooperative localization Ii29i to cooperative networks, provide 
insights into the network localization problem, and can guide 
the design and deployment of location-aware networks. 

The rest of the paper is organized as follows. Section |ll] 
presents the system model and the concept of SPEB. In Section 
HID we apply the notion of EFI to derive the SPEB. Then in 
Section |IV] we provide a geometric interpretation of EFIM 
for localization and derive scaling laws for the SPEB. Finally, 
numerical results are given in Section |V] and conclusions are 
drawn in the last section. 

Notation: The notation IEx{-} is the expectation operator 
with respect to the random vectors x; A ;^ B and A )^ B 
denote that the matrix A — B is positive definite and positive 
semi-definite, respectively; tr{ } denotes the trace of a square 
matrix; [-J^ denotes the transpose of its argument; [■]nxnk 
denotes the fcth nX7i submatrix that starts from element n{k — 
1 ) + 1 on the diagonal of its argument; [ • 1 ^ ^ ^ ^ denotes a 
submatrix composed of the rows ri to r2 and the columns ci 
to C2 of its argument; and || • || denotes the Euclidean norm of 
its argument. We also denote by /(x) the probability density 
function (PDF) /x(x) of the random vector X unless specified 
otherwise. 

II. System Model 

In this section, we describe the wideband channel model 
and formulate the localization problem. We briefly review the 
information inequaUty and the performance measure called 
SPEB. 

A. Signal Model 

Consider a synchronous network consisting of A^b anchors 
(or beacons) and A'a agents with fixed topology^ Anchors have 
perfect knowledge of their positions, while each agent attempts 
to estimate its position based on the waveforms received 
from neighboring nodes (see Fig. [U. Unlike conventional 
localization techniques, we consider a cooperative setting, 
where agents utilize waveforms received from neighboring 
agents in addition to those from anchors. The set of agents is 
denoted by = {1,2,..., N.^}, while the set of anchors is 
A^b = {iVa + 1,N^ + 2,...,N^ + Nb}. The position of node 
k is denoted by = [xk 2/fc]^0 Let (f)kj denote the angle 
from node k to node j, i.e., 

, . -1 Vk- Vj 
<pkj = tan — - — , 

•^k -^j 

'We consider synchronous networks for notional convenience. Our ap- 
proach is also vaUd for asynchron ous n etwo rks, where devices employ round- 
trip time-of-flight measurements 1251 . 1401 . 

''For convenience, we focus on two-dimensional localization where pj; 6 
R^, and we will later mention extensions to three-dimensional localization. 
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and Qfej ^ [cos(/)fcj sinc/jfcj Y denote the corresponding unit 
vector. 

The received waveform at the fcth agent (fc G N^) from the 
jth node (j e Wj, U 7Va\{fc}) can be written as jm, 11411 



(t) = E ~ )+^kjit). t e [ 0, To, 



(1) 



(=1 



where s{t) is a known wideband waveform with Fourier 
transform S{f), aj:j and are the amplitude and delay, 
respectively, of the Ith pathOLfej is the number of multipath 
components, Zkj (t) represents the observation noise, modeled 
as additive white Gaussian processes with two-sided power 
spectral density Nq/2, and [OjTob) is the observation interval. 
The relationship between the positions of nodes and the delays 
of the propagation paths is 



^(0 



1 



iPfe 



jeKuK\{k}, (2) 



where c is the propagation speed of the signal, and b^j^j > 
is a range bias induced by non-line-of-sight (NLOS) propaga- 
tion. Line-of-sight (LOS) signals occur when the direct path 
between nodes k and j is unobstructed, such that 6^^'' = 0. 



B. Error Bounds on Position Estimation 

We first introduce 6 as the vector of unknown parameters. 



e 



where P consists of all the agents' positions 



T T 
Pi P2 



T 



and dk is the vector of the multipath parameters associated 
with the waveforms received at the fcth agenj^ 



fc,fc+l 



''k,N.^+Nb 



in which Kkj is the vector of the multipath parameters asso- 
ciated with rfcj(t)n| 



'kj 



^kj 



Secondly, we introduce r as the vector representation of all 



the received waveforms, given by r 
where 



rT IT 



^k.N.^+Nb 



^We consider the general case where the wideband channel is not neces- 
sarily reciprocal. Our results can be easily specialized to the reciprocal case, 
where we have Li.,- = L,!., o!'? = and r,''.' = r'P hence bj'^ = b'^, 
for 1 = 1,2,.. .,Lkj. 

*In cases where the channel is reciprocal, only half of the multipath 
parameters are needed. Without loss of generality, we only use | a^'j , b^'j : 

k,j &Mi,k> j|. 



The bias 



for LOS signals. From the perspective of Bayesian 
estimation, it can be thought of as a random parameter with infinite a priori 
Fisher information 1291 . 



in which r^j is obtained from the Karhunen-Loeve (KL) 
expansion of rkj{t) B2l . 1431 . We tacitly assume that when 
nodes j and fc cannot communicate directly, the corresponding 
entry r^j is omitted in r. 

We can now introduce an estimator 6 of the unknown 
parameter 6 based on the observation r. The mean squared 
error (MSB) matrix of 6 satisfies the information inequality 

i-im 



E, 



(3) 



where Jg is the Fisher information matrix (FIM) for 00 given 
by 



ln/(r,0) 



(4) 



in which /(r, 6) is the joint PDF of the observation r and 
the parameter vector 0. For an estimate of the fcth agent's 
position, equation (|3]l implies that 



IEr,e {(Pfe - Pfc)(Pfc - Pfc)^} ^ [Je^] 



2x2, k ■ 



One natural measure for position accuracy is the average 
squared position error Er e {||pfc — pfc|p}, which can be 
bounded below by V{pk) defined in the following. 

Definition 1 (Squared Position Error Bound l[29]l ): The 
squared position error bound (SPEB) of the fcth agent is 
defined to be 



P(p,)^tr{[J-]^^^J 



Since the error of the position estimate pk — Pk is a vector, 
it may also be of interest to know the position error in a partic- 
ular direction. The directional position error along a given unit 
vector u is the position error projected on it, i.e., u^(pfe — p^), 
and its average squared error Er,e {||u^(Pa: ^ Pfc)lP} can be 
bounded below by Vipk] u) defined in the following]^ 

Definition 2 (Directional Position Error Bound): The 
directional position error bound (DPEB) of the fcth agent 
with constraint (p^ — p^) = is defined to be 

where u, u_l e are unit vectors such that (u, u^) = 0. 

Proposition 1: The SPEB of the fcth agent is the sum of 
the DPEBs in any two orthogonal directions, i.e.. 



r(pk) ^ V{pk;u) +V{pk;u^) ■ 



Proof: See Appendix lAl 



(5) 



□ 



*With a slight abuse of notation, E[..e{-} in (3j and will be used 
for deterministic, random, and hybrid cases, with the understanding that the 
expectation operation is not performed over the deterministic components of 
1431 , 1441 . Note also that for the deterministic components, the lower bound 
is valid for their unbiased estimates. 

'in higher dimensions, this notion can be extend to the position error in 
any subspaces, such as a hyperplane. 
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C. Joint PDF of Observations and Parameters 

Evaluation of (|4|i requires knowledge of the joint distri- 
bution f{r,9). We can write f{r,6) = f{r\9)f{9), where 
f{r\0) is the likelihood function, and f{6) is the a priori 
distribution of the parameter 00 In this section, we describe 
the structure of both functions in detail. 

Since the received waveforms rkj (t) are independent con- 
ditioned on the parameter 6, f{r\6) can be expressed as 1421 . 

ma 

/(rie)=n n /Ki^)' 



where 



firkj\0) (X exp 




rkj (t)^ ail s 



E 

1=1 



1=1 



dt 



dt 



(7) 

When the multipath parameters Kkj are independent condi- 
tioned on the nodes' positions^ f{0) can be expressed as 



fiO) 



/(P) 
/(P) 



(8) 



n 

n n /('^^•^■ip 

fceM jeAfbUMMk} 
where /(P) is the joint PDF of all the agents' positions, 
and f{Kkj\P) is the joint PDF of the multipath parameters 
Kkj conditioned on the agents' positions. Based on existing 
propagation models for wideband and UWB channels lfT4l . 
||251 . the joint PDF of the channel parameters can be further 
written as 



f{Kk,\P)^f{Kk,\dk,), (9) 

where dkj = \\pk - Pj\\ for fc e and j e A/i, U Afa\{k}. 
Combining dH) and (|6]l leads to 



ln/(r,0)= ^ [ln/(rfc,|0)+ln/K,|P) 



E 

keKje^fh 

J2 E [ln,/(rfc,|0)+ln/(Kfe,|P) 

keAf;]£K\{k} 



ln/(P) 



(10) 



where the first and second groups of summation account for 
the information from anchors and that from agents' cooper- 
ation, respectively, and the last term accounts for the infor- 
mation from the a priori knowledge of the agents' positions. 
This implies that the FIM for in (l4]i can be written as 
Je = + Je + where Jg, J^, and Jg correspond to 
the localization information from anchors, agents' cooperation, 
and a priori knowledge of the agents' positions, respectively. 



'"When a subset of the pai'ameters are deterministic, they are eliminated 
from f{e). 

"This is a common model for analyzing wideband communication, unless 
two nodes are close to each other so that the channels from a third node 
to them are coiTelated. Our analysis can also account for the con'elated 
channels, in which case the SPEB will be higher than that corresponding 
to the independent channels. 



III. Evaluation of FIM 

In this section, we briefly review the notion of EFI ll29l and 
apply it to derive the SPEB for each agent. We consider both 
the cases with and without a priori knowledge of the agents' 
positions. We also introduce the concept of RI, which turns 
out to be the basic building block for the EFIM. 

A. Equivalent Fisher Information Matrix and Ranging Infor- 
mation 

We saw in the previous section that the SPEB can be 
obtained by inverting the FIM Je in (|4|i. However, Je is a 
matrix of very high dimensions, while only a much smaller 
submatrix [Jq ^] 2Ar. x2Af. interest. To gain insights into 

localization problem, we will employ the notions of EFIM and 
RI 1291 . For the completeness of the paper, we briefly review 
the notions in the following. 

Definition 3 (Equivalent Fisher Information Matrix): 
Given a parameter vector 6 = [6l O^Y and the FIM Je of 
the form 



Je = 



A B 

BT C 



where 6 G R^, Oi e M", A e M"^", B e ^nx{N~n)^ 
and C € ^(N~n)x(N~n) ^jjjj „ < the equivalent Fisher 
information matrix (EFIM) for 6i is given by 



Je(0i) = A-BC-^B^ 



(11) 



Note that the EFIM retains all the necessary information to 
derive the information inequality for the parameter Oi, in a 
sense that [Jg^Jnxn = [Je(^i)]~^, so that the MSE matrix 
of the estimates of 6i is "bounded" below by [ Je(0i) ] ^. The 
right-hand side of (fTTT i is known as the Schur's complement 
of matrix A 145 L and it has been used for simplifying the 
CRBs ED, 132], ||46l. 

Definition 4 (Ranging Information): The ranging informa- 
tion (RI) is a 2 X 2 matrix of the form A Jr(0), where A is a 
nonnegative number called the ranging information intensity 
(RII) and the matrix Ji((/)) is called the ranging direction 
matrix (RDM) with the following structure: 



J,(</>) 



cos q> 
cos (h sin ( 



cos (p sm ( 



The RDM Ji (0) has exactly one non-zero eigenvalue equal 
to 1 with corresponding eigenvector q = [cos0 sinc/)]^, 
i.e., Ji(0) = qq^- Thus, the corresponding RI is "one- 
dimensional" along the direction 0. 

B. Fisher Information Analysis without A Priori Position 
Knowledge 

In this section, we consider the case in which a priori 
knowledge of the agents' positions is unavailable, i.e., /(P) 
is eliminated from ([8]l. We first prove a general theorem, 
describing the structure of the EFIM, followed by a special 
case, where there is no a priori knowledge regarding the 
channel parameters. 
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Theorem 1: When a priori knowledge of the agents' posi- 
tions is unavailable, and the channel parameters corresponding 
to different waveforms are mutually independent, the EFIM for 
the agents' positions is a 2N^ x 2Ns, matrix, structured as (fT2] i 
at the bottom of the page, where J^(pfc) and C^j can be 
expressed in terms of the RI: 

and 

Ckj = Cjk = (Afcj + Xjk) Ji-(0fcj) , 

with Xkj given by (l35T l in Appendix |B] 

Proof: See Appendix IbI □ 
Remark 1: We make the following remarks. 

• To obtain the SPEB of a specific agent, we can apply 
EFI analysis again and further reduce Je(P) into a 2 x 2 
EFIM. 

• The RI is the basic building block of the EFIM for 
localization, and each RI corresponds to an individual 
received waveform. The RJI Xkj is determined by the 
power and bandwidth of the received waveform, the 
multipath propagation, as well as the a priori channel 
knowledge. Note that each received waveform provides 
only one-dimensional information for localization along 
the angle 4>ky 

• The EFIM Je(P) can be decomposed into localization 
information from anchors and that from agents' cooper- 
ation. The former part is represented as a block-diagonal 
matrix whose non-zero elements are J^(pfe), for the kt\\ 
agent, and each J^(p/;) is a weighted sum of RDMs 
over anchors. Hence the localization information from 
anchors is not inter-related among agents. The latter 
part is a highly structured matrix consisting of RIs 
Ckj- Hence the localization information from agents' 
cooperation is highly inter-related. This is intuitive since 
the effectiveness of the localization information provided 
by a particular agent depends on its position error. 

Theorem 2: When a priori knowledge of the agents' posi- 
tions and the channel parameters is unavailable, the EFIM for 
the agents' positions is a 2N^ x 2N^ matrix, structured as in 
(fT2l i shown at the bottom of the page, with the RII Xkj given 
by 



A 



'8^2^Vc' • (1 - Xkj) SNR^^y , LOS signal, 
, NLOS signal, 



where j3 is the effective bandwidth of transmitted waveform 

s[t) 



13 = 



1/2 



SNR^^.^ is the SNR of the first path in rfej (t): 



No 



(13) 



and < Xkj < 1 is called the path-overlap coefficient, which 
depends on the first contiguous-cluste][3 in LOS signals. 

Proof: See Appendix O □ 
Remark 2: We make the following remarks. 

• The theorem shows that when a priori knowledge of 
channel parameters is unavailable, the NLOS signals do 
not contribute to localization accuracy, and hence these 
signals can be discarded. This agrees with the previous 
observations in JS), ||3T1 , ||32| although the authors con- 
sidered different models. 

• For LOS signals, the RII is determined by the first 
contiguous-cluster |[29l , implying that it is not necessary 
to process the latter multipath components. In particular, 
the RII is determined by the effective bandwidth /?, the 
first path's SNR, and the propagation effect characterized 

by Xk3- 

• Since Xkj > 0, path-overlap always deteriorates the 
accuracy unless Xkj = 0, in which the first signal compo 
nent s{t 



s{t 



for I > 1. 



t;„') does not overlap with later components 



C. Fisher Information Analysis with A Priori Position Knowl- 



We now consider the case in which the a priori knowledge 
of the agents' positions, characterized by /(P), is available. 
We first derive the EFIM, based on which we prove that 
agents and anchors can be treated in a unified way under this 
framework. We then present a special scenario in which the 
a priori knowledge of the agents' positions satisfies certain 
conditions so that we can gain insights into the EFIM. 

Theorem 3: When a priori knowledge of the agents' posi- 
tions is available, and the channel parameters corresponding 

'^The first contiguous-cluster is tlie first group of non-disjoint paths. Two 
patlis that arrive at time T j an d rj are called non-disjoint if | — | is less 
than the duration of s{t) 1291 . 



Je(P) 



J^(Pl)+ E Ci,, -Ci,2 

-Ci,2 J^(P2)+ E C2.J 

jeA4\{2} 



(12) 



-Cl 



^2,Af„ 



3<£K\{N,} 
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to different waveforms are mutually independent, the EFIM 
for the agents' positions is a 2Ns, x 2N.^ matrix, given b\F^ 



Je(P)=J^(P)+Je^(P)+Hp, 



(14) 



where 



2fc,2m-l:2m \q k^m, 



[ Je (-P) ]2fc-l:2fc,2m- 



1:2m 

+ Rfc(r™fc)] , 



k = m , 
k m , 



and 



Hp =Ei 



52 



in/(p; 



with Rfe(rfcj) e M^^^ given by ( fTSl ) shown at the bottom of 
the page. Block matrix (x, y) in ( fTSl l is defined as (l27T i 
in Appendix IbI 

Proof: See Appendix ID] □ 

Remark 3: The EFIM for agents' positions is derived in 
( fT4] i for the case when a priori knowledge of the agents' 
positions is available. Compared to ( fT2] i in the Theorem [T] 
the EFIM in ( fT4l i retains the same structure of the localization 
information from both anchors and cooperation, except that all 
RIs in Theorem[3]are obtained by averaging the 2x2 matrices 
over the possible agents' positions. In addition, the localization 
information from the position knowledge is characterized in 
terms of an additive component Hp. This knowledge improves 
localization because Hp is positive semi-definite. 

Based on the result of Theorem[3] we can now treat anchors 
and agents in a unified way, as will be shown in the following 
theorem. 

Theorem 4: Anchors are equivalent to agents with infinite 
a priori position knowledge in the following sense: when the 
kth agent has infinite a priori position knowledge, i.e., Hp^ = 
limj2_^g^ diag t^}, then 



Je(Ps) = [Je(P)] 



k ' 



"Note that Jc(P) in 03) does not depend on any particular value of the 
random vector P, whereas Jc(P) in )12t is a function of the deterministic 
vector P. 



where P;, is the vector P without rows 2fc — 1 to 2k, and 
[ Je(P) ]j. is the matrix Je(P) without rows 2fc — 1 to 2k and 
columns 2k — 1 to 2k. 

Proof: See Appendix |E] □ 

Remark 4: The theorem shows mathematically that agents 
are equivalent to anchors if they have infinite a priori position 
knowledge, which agrees with our intuition. As such, it is 
not necessary to distinguish between agents and anchors. This 
view will facilitate the analysis of location-aware networks and 
the design of localization algorithms: every agent can treat the 
information coming from anchors and other cooperating agents 
in a unified way. 

The general expression of the EFIM for the case with a 
priori position knowledge is given in (fT4l) . which is much more 
involved than that for the case without position knowledge in 
(fT2] |. However, in the special case when 



Ep {5(P)} = g(Ep{P}), 



(16) 



for the functions q(-) involved in the derivation of the EFIM 



4^' 



(see Appendix |D]io we can gain insight into the structure of 
the EFIM as shown by the following corollary. 

Corollary 1: When the a priori distribution of the agents' 
positions satisfies (fTSl i. and the channel parameters corre- 
sponding to different waveforms are mutually independent, 
the EFIM for the agents' positions is a 2A^, x 2iVa matrix, 
structured as ( fTTI l shown at the bottom of the page, where 
J^(Pfc) and Ckj can be expressed in terms of the RJ: 

J^(P/c) = X! ^fej Jr(0fej), 



and 



Cfcj = Cjk = (Afej + \jk) 3y{(t)kj) . 



where P = Ep {P}, Xkj is the RII given in i35[ evaluated at 
P, and (f>kj is the angle from to pj. 

Proof: See Appendix ID] □ 



D. Discussions 

We will now discuss the results derived in the previous 
sections. Our discussion includes 1) the EFIM for the agents 

This occurs when every agent's a priori position distribution is concen- 
trated in a small area relative to the distance between the agent and the other 
nodes, so that g{P) is flat in that area. 



R.fc(rfcj) = Ep {$fej {dkj,dkj) qkj qlj} - Ep {qkj *fej(4j,Pfc)} Ep {^kj{K,kj , K,kj)} ^ Ep {^kj{Pk, dkj) qlj} (15) 



Je(P) 



J^(Pi) + E Ci,,- 

jeK\{i} 



-Ci, 



J^(P2)+ E C2,, 
J6A4\{2} 



^2,Ar. 



Cp 



(17) 
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in non-cooperative localization, 2) an application of the co- 
operative localization to tracking, 3) a recursive method to 
construct an EFIM for large networks, and 4) the extension to 
three-dimensional scenarios. 

1) Non-Cooperative Localization: When the agents do not 
cooperate, the matrices corresponding to the agents' coopera- 
tion in (fTZt in Theorem[T]and ( fTTI i in Corollary [T] are discarded. 
In particular, the EFIM Je(P) in Theorem [T] reverts to 

Je(P)=diag{j^(pi),J^(p2), ... ,J^(pjvj}, 

and hence the 2x2 EFIM for the fcth agent is equal to 
Je(Pfe) Jf'(Pfc)- Similarly, the EFIM Je(P) in Corollary 
[T] reverts to 

Je(P)=diag{jf(pi),J^(p2), ... ,jA(p^j}+Hp. 

Furthermore, when the agents' positions are independent a 
priori. Hp = diag {Hpi , Hp^ , . ■ . , 3p„^ } and the 2x2 EFIM 
for the fcth agent can be written as Je(pfc) = J^(pfc) + Sp^.. 

2) Spatial vs. Temporal Cooperation for Localization: 
Rather than multiple agents in cooperation, a single agent can 
"cooperate" with itself over time. Such temporal cooperative 
localization can easily be analyzed within our framework, as 
follows. 

Consider a single agent moving in sequence to N different 
positions according to piecewise linear walk and receiving 
waveforms from neighboring anchors at each position. The 
N positions can be written as P = [ p j pj • • • ] , and 
we can consider the scenario as N agents in cooperation. The 
likelihood of the observation is 

/ (r,di©) = n n fi^kj\Pk,p,) n /(4iPfc,Pfc+i) , 



k=l 



where d ~ [di d2 ■■■ dAr-i]^ in which dk is the 
measurement of the distance dk = \\pk — Pfe+i|| between p^ 
and Pfc+iIH By applying Theorem [5 we have the EFIM for 
P as Je(P) = Jf'(P) + J^(P) where 

J^(P) = diag { J^(pi), J^(P2), . . . , J^(Piv)} , 

'^We assume that the agent has other navigation devices, such as inertial 
measurement unit (IMU), odometer, or pedometer, to measure the distance 
between positions. 



and Jg (P) is given by (fTsT l shown at the bottom of the page, 
in which C^^fe+i = i^k ^rifpkM+i) with (j>k,k+i denoting the 
angle from p^ to Pfe+i and 

( f - 

Vk =IEd f \ dk\pk,Pk+i 

By further applying the notion of EFI, we can obtain the EFIM 
Je(pfc) for each position p^. Note that this analysis can be 
extended to cooperation among multiple mobile agents over 
time, so that both cooperation over space and time are explored 
simultaneously. 

3) Recursive Formula for EFIM: The structure of the EFIM 
in ( fT2] i and ( fTTb enables us to extend the EFIM when agents 
join or leave the cooperative network. We will develop a 
recursive formula to construct the EFIM in the following. 

Consider a network with n agents in cooperation without a 
priori knowledge of their positions, and the EFIM for agents' 



positions Je(P,i) where P„ = [p| • • • p\Y can be obtained 



by ( fT2b . If a new agent enters the cooperative network, then 
the EFIM for the n + 1 agents is given by ( fT9b . shown at 
the bottom of the page, where Ja^m+i is the EFIM for the 
(n + l)th agent corresponding to the localization information 
from anchors, M„ „+i is the localization information from 
the cooperation between the (n + l)th agent and the other n 
agents, given by 

M„ = diag {Ci , C2,n+1 J • ■ • J C„_„+i} , 

and K„ e ]g2nx2 given by 



Kn — I2x2 I2x2 ••■ I2x2 



Note that when the a priori knowledge of the agents' positions 
is available, we need to consider the contribution of Hp, and 
the EFIM for the n + 1 agents can be constructed in a similar 
way. 

Similarly, when a certain agent, say fc, leaves the network, 
we need to eliminate rows 2fc — 1 to 2k and columns 2fc — 1 
to 2fc in Je(P„), as well as subtract all corresponding Ckj for 
j e Mi\{k} from the diagonal of Je(P„). 

4) Extension to 3D Localization: All the results obtained 
thus far can be easily extended to the three-dimensional 
scenario, in which pfc = [xk Vk ZkY ■ The SPEB of the 
fcth agent is defined as Vipk) — [Ja^Iq , i.- Following the 
steps leading to (fl2l) and dlTl ). we can obtain a corresponding 



Cl,2 



Je (P) = 



Cl,2 + C2,3 
— C2,3 



■^2,3 



(18) 



— CjV-l,JV Cat-i jv 



Je(P, 



n+lj 



Je(P„) + M„,„+i -M„,„+i K„ 



(19) 



8 



IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. X, NO. Y, MONTH 2010 





X* 

X 

^ ^""-^ 






/ \ 

/ \ 


.-"^^ / 


/ y 
1 / 


\ / 


1 / 
\ / 


\ 






y ^ 


\ 

\ 



Fig. 2. Geometric inteipretation of the EFIM as an information ellipse. In 
the rotated coordinate system (rotated over an angle the major and minor 
axes of the ellipse are given by yfjl and ^Jt], respectively. 




X 



Fig. 3. Updating of the information ellipse for non-cooperative localization. 
The original information ellipse of the agent is characterized by F(/j, 0). 
The RI from an additional anchor is given by F{i^, 0, (/>). The new information 
ellipse of the agent then grows along the direction </>', but not along the 
orthogonal direction. The new information ellipse corresponds to Y(Ji,f\,'d). 



3iVa X 3iVa EFIM involving the RDMs Jj. ((/j^^ , ) for fc e 
and i e A/b U A/'a, where 

with (ys and denoting the angles in the spherical coordinates, 

and q = [ cos COS (/) sincpcos(/) sint/)]^. 



IV. Geometric Interpretation of EFIM for 
Localization 

In this section, we present a geometric interpretation of 
the EFIM for localization. This interpretation not only pro- 
vides insights into the essence of localization problems, but 
also facilitates the analysis of localization systems, design 
of localization algorithms, and deployment of location-aware 
networks. We begin with the non-cooperative case, and then 
extend to the cooperative case. Based on these results, we 
derive scaUng laws of the SPEB for both non-cooperative and 
cooperative location-aware networks. 

A. Interpretation for Non-Cooperative Localization 

When an agent only communicates with neighboring an- 
chors, the EFIM can be written a 



Je(p)= AjJr(<^,)=U^ 



/i 
?/ 



u 



■i5 1 



(20) 



where /i and rj are the eigenvalues of Je(p), with n > rj, and 
is a rotation matrix with angle d, given by 



cos d — sin 
sin i} cos d 



The first and second columns of are the eigenvectors 
corresponding to eigenvalues fi and 77, respectively. By the 
properties of eigenvalues, we have 

^ + 77 = tr{Je(p)} =J2^i- 

'*To simplify the notation, we will suppress the agent's index in the 
subscript. 



Note in ( |20l i that Je(p) depends only on /^i, 77, and i), and we 
will denote Je(p) by F(/i, 77,7?) when needed. 

Proposition 2: The SPEB is independent of the coordinate 
system. 

Proof: See Appendix IF] □ 
Remark 5: The proposition implies that if we rotate the 
original coordinate system by an angle -d prescribed by (l20b 
and denote the agent's position in the new coordinate by p*, 
then the SPEB is 



^(P) = Vip*) = tr 



H 
?7 



1 
/i 



The EFIM in the new coordinate system is diagonal, and thus 
the localization information in these new axes is decoupled. 
Consequently, the SPEB is also decoupled in these two or- 
thogonal directions. 

Definition 5 (Information Ellipse): Let J be a 2 x 2 positive 
definite matrix. The information ellipse of J is defined as the 
sets of points x e such that 



X J-^x^ 



1. 



Geometrically, the EFIM in ( [20] i corresponds to an informa- 
tion ellipse with major and minor axes equal to y/JI and y/fj, 
respectively, and a rotation -d from the reference coordinate, as 
depicted in Fig.|2] Hence, the information ellipse is completely 
characterized by fi, 77, and 7?. Note that the RI is expressed 
as XJi-{(f>) = F(A, 0,(/>), and it corresponds to a degenerate 
ellipse. In the following proposition, we will show how an 
anchor contributes to the information ellipse of an agent. 

Proposition 3: Let Je(p) = F(/.i, 77,7?) and V{p) denote 
the EFIM and the SPEB of an agent, respectively. When that 
agent obtains RJ F(i/, 0, </>) from a new anchor, the new EFIM 
for the agent will be 



Je(p)=F(/2,77,^) 
= F(a7,77,79) 



F(i.,0,, 
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where the parameters for the new information ellipse are 

^ ^ ^ ^ + i - r; + cos 20'] V 20' 



and 



2 2 

+ + ly 1 

2 2 
1 



\J ri + u COS 20']^ + sin^ 20' 

^ sin 20' 

V = V A — arctan , 

2 ji — 1] + V cos 20' 

with 0' = — I?. Correspondingly, the new SPEB becomes 



V{V) = \ + \ = 



11 + 11+ u 



A ry IJirj + V [rj + {^1 - vi) svcl 0'J 



(21) 



Remark 6: The geometric interpretation for the proposition 
is depicted in Fig. [3] For a fixed Rll ly, we see from ( ITTT i that 
P(p) can be minimized through 0' (equivalently, through 0) 
in the denominator: 



min V{p) = 



and the minimum is achieved when = ± 7r/2. In such 
a case, the anchor is along the direction of the eigenvector 
corresponding to the smallest eigenvalue rj. Observe also that 
the denominator in (l2Tl i is equal to /i • fj, which is proportional 
to the squared area of the new information ellipse correspond- 
ing to Je(p). Hence, for a fixed i/, the minimum SPEB is 
achieved when the new anchor is along the minor axis of 
the information ellipse corresponding to Je(p). Equivalently, 
this choice of anchor position maximizes the area of the new 
information ellipse. 

On the other hand, the maximum SPEB occurs when the 
anchor is along the direction of the eigenvector correspond- 
ing to the largest eigenvalue fi, i.e., the major axis of the 
information ellipse corresponding to Je(p). Equivalently, this 
minimizes the area of the new information ellipse, and thus 

max V[l>) = — — ■ — - , 
<t> ri[fi + iy) 

and the maximum is achieved when = ± tt. Note also that 

1 

M 

min^ V{p), and the 



< V{p) < Vip) , 



where the left-hand side l/fJ. ^ lim^^ 
right-hand side 7'(p) = limi^^o^(p)- 



B. Interpretation for Cooperative Localization 

The EFIM for all the agents in cooperative location-aware 
network is given respectively by iT% and (fT2l i for the cases 
with and without a priori position knowledge. Further applying 
the notion of EFl, one can obtain the EFIM for individual 
agents. In general, the exact EFIM expression for the individ- 
ual agents is complicated. However, we can find lower and 
upper bounds on the individual EFIM to gain some insights 
into the localization problem. 

Proposition 4: Let J^(pfc) = F(/ifc, 77^, i^^) denote the 
EFIM for agent k that corresponds to the localization infor- 
mation from anchors, and let Ckj = F{i'kj, 0, (j>kj) denote the 



Agent 2 



J^{P2) 




Agent 1 







1\ 


)''- - 



Je(Pl) 



'J^(Pl) 



Fig. 4. Updating of the information ellipse for cooperative localization. 
Based on the anchors, the feth agent has information J^(pj.). The cooperative 
information between the two agents is given by Ci 2 = F{u,0, (j>i^2)- The 
total EFIM for agent 1 is then Je{pi) = J*{pi) + Ci,2Ci,2- The new 
information ellipse grows along the line connecting the two agents. 



RI for that agent obtained from cooperation with agent j. The 
EFIM Je(Pfc) for agent k can be bounded as follows: 



where 



Je (Pfc) ^ Je(Pfc) ^ Jl^iPk) , 



J^(Pfc) = Je^(Pfc) + E ^k.Ckj, (22) 



Je"(Pfe)=J^(Pfc)+ E 



C 



kj 



(23) 



with coefficients < ^Jr^ < < 1 given by (gUi and ( |46] |. 
Proof: See Appendix |F] □ 

Remark 7: The bounds for the EFIM can be written as 
weighted sums of RIs from the neighboring nodes, and such 
linear forms can facilitate analysis and design of location- 
aware networks. Moreover, it turns out that ^j:^ = when 
there are only two agents in cooperation, leading to the 
following corollary. 

Corollary 2: Let J^(pi) = F(^i, 771, 7?i) and J^(p2) = 
Y{p.2,r]2,'&2) denote the EFlMs for agent 1 and 2 from 
anchors, respectively, and let C1.2 = F(i/i 2, 0, 0i_2) denote 
the RI from their cooperation. The EFlMs for the two agents 
are given, respectively, by (see also Fig. |4|i 



and 



where 



and 



Je(pi) = J^(Pl) + 6,2 J^l,2 Jr(01,2) 



Je(p2) = J^(P2) + 6,1 ^1,2 Jr(01,2) 



6,2 = 



1 



6,1 = ^ 



1 + Vi^2 A2(01,2) ' 
1 



Via Ai 
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with 

Afc(0i,2) = qL [J^(Pfc)]"'qi2, 

for fc = 1,2. 

Remark 8: The results follow directly from Proposition |4] 
We make the following remarks. 

• Cooperation provides agent 1 with RI ^1,2 1^1,2 Jr(</'i,2) 
with < ^1,2 < 1. Hence agent 1 obtains a RII 
■^1,2 1^1,2 from cooperation instead of the full RII j/i,2. 
This degradation in RII is due to the inherent uncertainty 
of the second agent's position. We introduce the effective 

RII i>i_2 = 6,2 1^1,2- 

• The effective RII has the following geometric interpreta- 
tion. The value A2(0i.2) is the DPEB of agent 2 (based 
solely on the anchors) along the angle 0i 2 between the 
two agents. This implies that the larger the uncertainty 
of agent 2 along the angle 0i,2, the less effective coop- 
eration is. For a given A2{4>i.2), the effective RII i>i^2 
increases monotonically with i^i 2, and has the following 
asymptotic Umits: 

lim Di^2 = , 

lim !>i.2 1/A2(01,2) • 

Hence the maximum effective RII that agent 2 can 
provide to agent 1 equals the inverse of the DPEB of 
agent 2 (based solely on the anchors) along the angle 
cf>i^2 between the two agents. 

• When i) the two agents happen to be oriented such that 
(pi ^2 = 1^2, and ii) agent 2 is certain about its position 
along that angle (/i2 = +00), then A2((/)i,2) = and 
Je(pi) = J^(Pi) + Ci,2, i.e., agent 2 can be thought 
of as an anchor from the standpoint of providing RI 
to agent 1. From this perspective, anchors and agents 
are equivalent for localization, where anchors are special 
agents with zero SPEB, or equivalently, infinite J^(pfe) 
in all directions. 

C. Scaling Laws for Location-Aware Networks 

In this section, we derive scaling laws of the SPEB for both 
non-cooperative and cooperative location-aware networks. 
Scaling laws give us insight into the benefit of cooperation 
for localization in large networks. As we will see, agents and 
anchors contribute equally to the scaling laws for cooperative 
location-aware networks. 

We focus on two types of random networks: dense networks 
and extended networks 1471 . Il48l . In both types of networks, 
we consider the iVb anchors and agents randomly located 
(uniformly distributed) in the plane. In dense networks, adding 
nodes increases the node density, while the area remains 
constant. In extended networks, the area increases proportional 
to the number of nodes, while both the anchor and the 
agent densities remain constant. Without loss of generality, 
we consider one round of transmission from each node to 
another. All transmission powers are the same, while large- and 
small-scale fading can be arbitrary. Medium access control is 
assumed so that these signals do not interfere with one another 



Definition 6 (Scaling of SPEB): Consider a network with ?i 
nodes randomly located in a given area. We say that the SPEB 
of individual agents scales as 0(/(n)) for some function 
/(n), denoted by ^(p) £ 8(/(n)), if there are deterministic 
constants < ci < C2 < +00 such that 

P {ci/(n) < P(p) < C2/(n)} = 1 - e(n) , (24) 

where lim„^oo ~ 0. 

Theorem 5: In dense networks, the SPEB of each agent 
scales as 0(l/iVb) for non-cooperative localization, and as 
0(l/(iVb + iVa)) for cooperative localization. 

Proof: See Appendix iGl □ 

Theorem 6: In extended networks with an amplitude loss 
exponent 60 the SPEB of each agent scales as 

re(l/logAfb), b=l, 
V{p)ele{i), b>i, 

[ 9(1^-1), 0<6<1, 

for non-cooperative localization, and 

re(l/l0g(7Vb + iVa)), fe-1, 

P(p)eie(i), fe>l, 

for cooperative localization. 

Proof: See Appendix iGl □ 
Remark 9: We make the following remarks. 

• In dense networks, the SPEB scales inversely proportional 
to the number of anchors for non-cooperative localization, 
and inversely proportional to the number of nodes for 
cooperative localization. The gain from cooperation is 
given by 8(1 + N^/N^), and hence the benefit is most 
pronounced when the number of anchors is limited. 
Moreover, it is proven in Appendix |G] that e{n) decreases 
exponentially with the number of nodes. 

• In extended networks with an amplitude loss exponent 
equal to 1, the SPEB scales inversely proportional to the 
logarithm of the number of anchors for non-cooperative 
localization, and inversely proportional to the logarithm 
of the number of nodes for cooperative localization. This 
implies that the SPEB in extended networks decreases 
much more slowly than that in dense networks, and the 
gain from cooperation is now reduced to 9(log(A'b + 
A^a)/ log iVb). Moreover, it is shown in Appendix IgI that 
e{n) decreases as cxp(— (log 71)^/8)/ log n. 

> In extended networks with an amplitude loss exponent 
greater than 1, the SPEB converges to a strict positive 
value as the network grows. This agrees with our intu- 
ition that as more nodes are added, the benefit of the 
additional nodes diminishes due to the rapidly decaying 
RII provided by those nodes. It can be shown that the 
SPEB converges to a smaller value in the cooperative 
case than that in the non-cooperative case, i.e., a constant 
gain can be obtained by cooperation. 

"Note that the amplitude loss exponent is b, while the con'esponding power 
loss exponent is 2b. The amplitude loss exponent b is environment-dependent 
and can range from approximately 0.8 (e.g., hallways inside buildings) to 4 
(e.g., dense urban environments) 1491 . 
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Fig. 5. Effective RII ^12 !^i,2 as a function of the RII 1^1,2, for Ja,2 = 
F(/i2 = 2,r)2 = 1,i?2 = 0), and different angle of arrival tj>i^2- 



Fig. 6. Typical network deployment of two sets of anchors (set I: squares, 
set II: diamonds) and A^a = 15 agents. The agents are distributed uniformly 
over the [—10, 10] X [—10, 10] map, while the locations of the anchors are 
controlled by D. 



V. Numerical Results 

In this section, we examine several numerical examples 
pertaining to cooperative localization and illustrate practical 
applications of our analytical results. 

A. Effective Ranging Information 

We first investigate the behavior of the effective RII 1)1,2 
from Corollary |2] when two agents cooperate. The effective 
RII 1/12 is plotted in Fig. |5]as a function of the RII 1^1 2 for 
J^(P2) = F(/i2 = 2,772 = 1,1^2 = 0) and various values of 
(t>i,2- The corresponding asymptotic limits are also plotted for 
large values of i'i^2- We observe that effective RII increases 
from to l/A2((/>i,2) as the RII z^i 2 increases. For a fixed 
RII, the second agent will provide the maximum effective 
RII at 01,2 = 7?2, along which angle the second agent has 
the minimum DPEB (i.e., l/fJ,2 = 0.5). On the other hand, 
the second agent will provide the minimum effective RII at 

2 = 192 ± 7''/2, along which angle the second agent has the 
maximum DPEB (i.e., 1/772 = !)■ 
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Fig. 7. The average SPEB as a function of the number of agents in the 
network for various anchor configurations (D = 10). 



B. Benefit of Cooperation 

We now consider the SPEB performance as a function 
of the number of agents for cooperative localization. The 
network configuration is shown in Fig.|6] The agents randomly 
(uniformly distributed) reside in a 20 m by 20 m area. There are 
two sets of anchors (shown as squares (set I) and diamonds (set 
II) in Fig.|6]l, with a configuration determined by the parameter 
D. Since fading does not affect the scaling behavior as shown 
Section IIV-CI we consider a network with signals that obey 
the free-space path-loss model for simplicity, so that the RII 

Figure [7] shows the average SPEB over all the agents as a 
function of the number of agents, obtained by Monte Carlo 
simulation, for D = 10. We see that as the number of agents 
increases, the average SPEB decreases significantly, roughly 
proportional to the number of agents. Note that the anchor 



configuration set II yields a lower SPEB. Intuitively, this is 
due to the fact that the anchors in set II (distance D from the 
center) cover the area better than the anchors in set I (distance 
\/2D from the center). 

Define the upper and lower approximations of agent fc's 
SPEB as 

7'"(p.)^tr{[jL(p,)]-^} 

and 

7'^(p.)^tr{[J>,)]-^}, 

where Jg(pfe) and J^(pa;) are given by ( |22] | and ( l23T l. re- 
spectively, in Theorem |4] Figure [8] shows the average ratio of 
the lower and upper approximations of the SPEB, obtained by 
Monte Carlo simulation, for anchor set I, set II, and both sets. 
When there are only two agents in cooperation, the bounds 
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Fig. 8. Ratio of upper and lower approximations of the SPEB, 'P^{p) and 
7'^(p), as a function of the number of agents for anchor set I, set II, and 
both. 



coincide, as we expect from Corollary [2] As the number of 
agents increases, the ratio deviates from 1, or equivalently 
the approximations become looser, due to the fact that upper 
approximation ignores more cooperative information, and the 
lower approximation considers more agents to be equivalent 
to anchors. Nevertheless, the ratio converges to a positive 
constant, implying that the upper and lower approximation 
decrease at the same rate in an asymptotical regime, as shown 
in the proof of Theorem |5] 

C. Anchor Deployment 

Finally, we investigate the effect of anchor deployment in 
more detail. We consider a scenario with = 15 agents. 
The anchor placement is controlled through D (see Fig. |6]l. 
Figure |9] shows the average SPEB as a function of D for 
different anchor configurations (set I, set II, and both sets). 
We see that the SPEB first decreases, and then increases, 
as a function of D. When D is close to 0, all the anchors 
are located closely in the middle of the area, and hence the 
RIs from those anchors to a particular agent are nearly in 
the same direction. This will greatly increase the error of 
each agent's position since every J^(pfc) is close to singular, 
resulting in poor overall SPEB performance. As the anchors 
begin to move away from the center, they provide RIs along 
different directions to each agent, which lowers the average 
SPEB. Then, as the distances of the anchors to the center 
increase further, the anchors become far away from more and 
more agents. Hence the RJI decreases due to the path-loss 
phenomena, and this leads to the increase in the average SPEB. 
Observe also that anchor set I is better than anchor set II for 
D < 7m. This is because, for a fixed _D < 7m, anchor set 
I can cover a larger area. For _D > 7 m, anchor set I suffers 
more from path-loss than anchor set II. 

For the sake of comparison, we have also included the 
average SPEB when 8 anchors are deployed 1) according to set 
I and II simultaneously, and 2) randomly on a [— 10 m, 10 m] x 




2 4 6 8 10 12 14 16 



D (m) 

Fig. 9. The mean SPEB with respective to anchor deployment. There are 
A^a = 15 agents. 

[—10 m, 10 m] area. The figure shows that intelligent anchor 
deployment can be beneficial compared to random deploy- 
ment, indicating the need for anchor deployment strategies. 

VI. Conclusion 

In this paper, we have investigated the fundamental limits on 
the localization accuracy for wideband cooperative location- 
aware networks. We have derived the squared position error 
bound (SPEB) by applying the notion of equivalent Fisher 
information (EFI) to characterize the localization accuracy. 
Since our analysis exploits the received waveforms rather 
than specific signal metrics, the SPEB incorporates all the 
localization information inherent in the received waveforms. 
Our methodology unifies the localization information from 
anchors and that from cooperation among agents in a canonical 
form, viz. ranging information (RI), and the total localization 
information is a sum of these individual RIs. We have put 
forth a geometrical interpretation of the EFIM based on eigen- 
decomposition, and this interpretation has facilitated the theo- 
retical analysis of the localization information for cooperative 
networks. We have also derived scaling laws for the SPEB 
in both dense and extended networks, showing the benefit 
of cooperation in an asymptotic regime. Our results provide 
fundamental new insights into the essence of the localization 
problem, and can be used as guidelines for localization system 
design as well as benchmarks for cooperative location-aware 
networks. 

Appendix A 
Proof of Proposition[T] 

Proof: The right-hand side of (|5]l can be written as 
7'(pfe;u)+7'(pfe;ui) 

= tr{uT [J,i]^^^^^u}+tr{ul [J,i]2,2,fe 
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where we have used the fact u + u i uT =1. 



□ 



Appendix B 
Proof of Theorem[T] 

We proceed in two steps: we first show that the EFIM is 
structured as in (fT2l i. and then derive the details of the RI. 

A. Derivation of the EFIM Structure 

When a priori knowledge of the agents' positions is unavail- 
able, the log-likelihood function in ( fTOt becomes 

ln/(r,K|P)=^ X] ln/(rfcj|Pfc,Pj,'«fci) 

feeM, jeA4uM,\{fe} 

+ ln/(Kfej|pfe,Pj) , (25) 

where k, denotes the vector of the channel parameters contain- 
ing all Kkj with k G and j e A/i, U A/'a\{fc}- For notational 
convenience, we now introduce 

(92ln/(r,M|P) 



*(x,y) ^Er 



92 



ln/(rfcj|pfe,Pj,Kfej) 
-Mn/(/«fcj|pfc,Pj) 



(26) 



(27) 



as well as 

Y(x, y, z) ^ *(x, y) [*(y, y)] #(y, z) , 

Tfcj(x,y,z) = *fcj(x,y) [*fcj(y,y)]"^ *fcj(y,z) • 

Since *(6ifc,0j) = for 7^ j, the EFIM for P can be 
derived as 

Je(P) = *(P,P)- ^ T(p,0fc,p) . (28) 

Structure 0/$ (P, P); Due to the structure in ( |25] |. we can 
express $(P,P) as 

* (p, p) = ^ 5] (p, p) + 5] p) 

fceXjGAAi, fceA/;,ieA4\{fc} 
= Ka + Kc , 

where Ka G M2Afax2JVj ^ block-diagonal matrix, consisting 
of 2 X 2 block matrices, given by 

/ EjeA4 (Pfc' Pfc) ' = : 



Structure 0/ Y ^P, fl'fe, P^ ; Since ^{Kki,Kkj) = for 
i 7^ J, we find that 



T P,6»fc,P 



5] 5] Y,,(P,K,„P) 

keM,jeMMk} 
^ Ma + Mc , 

where Ma G ]R2Afi,x2Af„ ^ block-diagonal matrix, consisting 
of 2 X 2 block matrices, given by 



[Ma ]2fe_i:2fc,2m-l:2m 



0. 



m , 
k ^ m . 



On the other hand, Mc G M^^"^^^" is also a block-matrix, 
consisting of 2 x 2 block matrices, given by (|30] | shown at the 
bottom of the page. 

Structure 0/ Je(P)." Combining these results, we find that 
the EFIM in (|28]| can be written as 



Je(P) = {Ka - Ma} + {Kc - Mc} , (31) 



from anchors 



from cooperation 



[K 



A J2/c-l:2fe.2m-l:2?i 



0. 



k ^ m . 

On the other hand, Kc G R2JVax2JV, ^j^^ ^ block-matrix, 
consisting of 2 x 2 block matrices, given by (|29] l shown at the 
bottom of the page. 



from which we obtain fT% . In ( fT2l i. 3^{pk) = 
J2jGN-k'^k{rkj) and Ckj = Cjfc = Rk{rkj) + Rki^jk) in 
which we have introduced the RI: 

Rfe(rA;j) = *fej (pfc,Pfc) - Tfcj (pfc,Kfcj,Pfe) ■ (32) 
Note that in the derivation, we used 

*fem (Pfe, Pm) = -*fcm (Pfc, Pfc) , 

and 

Yfem (PA;, '^fem, Pm) = km (P/c, '^fem, Pfe) • 

Since Je(P) in ( fT2] i can be expressed in terms of the RIs 
R-fe(rfe-j), for fc G A/'a and j G A/i, U A/'aXjfc}, we will examine 
next the details of the RIs. 

B. Details of the Ranging Information 

We now consider the detailed expression of the RI Rfc {vkj) 
in ( |32] ). We first introduce 



Sfej (x, y) = <^ - 



92 l n/(Kfcj-|Pfe,Pi) 
9x9yT 



and 



C J2fe-l;2fc,2m-l:2m, 



*fem (Pfe, Pm) + *mfe (Pfe, Pm) , k ^ m . 



(29) 



[M, 



C J2fc-l:2fe,2m-l:2n 



EjeAr,\{fe} [^/^i (Pfe'^^'^i'P'^) + ^J"*: (Pfe''*jfc'P'=)] ' = 

(pfc ) , k ^ m . 



(30) 
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a2ln/(rfcj|pfc,pj,/tfcj) 



(33) 



where kkj — [r^ 



(1) ^(2) ~(2) 



with a 



/c. 



a 



a 



kj 



From ^ and we note that dfcj = ||pfc — Pj|| and that 
/(rfcj|pfc,Pj , /«fcj) and / (/«fcj|pfc, Pj) only depend on Pfc,Pj 
through dkj- Using the chain rule, we have 



*/cj (Pfc,Pfc) 



^kj [dkj^dkj) 



and 



9pfc 

ddkj 
dpk 



dpi 



dd. 



k] 



Tkj {dkj,Kkj,dkj, 



and hence {vkj ) can be expressed as 

Rfe (r/cj) = ^k] {dkj,dkj) qkj qlj 

- Tfej (dfcj , fifcj , dkj ) Qfcj 
= Afej Qfej qL , 



(34) 



where q^j = ddkj/dpk = -ddkj/dpj = 
[cos (j)kj sin^fcj]^, and A^j is given by (l35T l shown at the 

bottom of the page, where Ikj = 



1 



1 



Appendix C 
Proof of TheoremE] 

Proof: When a priori channel knowledge is unavailable, 
we have Ekj{dkj,dkj) = 0, and 3^^ (4^ , k^j ) = 0. For 
NLOS signals, the RII in (l35l l becomes Xkj = since 
Hfcj {Kkj , Kkj ) = 0. For LOS signals, however, after some 
algebra, the RII becomes ( |36] | shown at the bottom of the 
page, where akj{K.kj, Kkj) = limj2^o^ diag O} since the 
Fisher information for known fej.^' = is infinity. To simplify 



we partition VP^j as 



kj 



where u^^ = Stt^/?^ SNR^^ obtained from ( |33] | through some 
algebra. As — > oo in (|36] |, we have 



^(1-X..)SNR« 



A. 



kj 



where 



Xkj 



87r2^2SNRa) 



(37) 



is called path-overlap coefficient 

We next show that only the first contiguous-cluster contains 
information for localization. Let us focus on Xkj ■ If the length 
of the first contiguous-cluster in the received waveform is Lkj, 
where 1 < Lkj < Lkj, we have 



''kj 



0' 



and *i. 



kj 











where kkj G M^l^,-:^ ^ ^{2L,^-i)x{2L^,-i) ^ ^j^^j ^ 
a block matrix that is irrelevant to the rest of the derivation. 
Hence ( [37] i becomes 



Xkj 



87r2/32SNR^^ 



which depends only on the first Lkj paths, implying that 
only the first contiguous-cluster of LOS signals contains 
information for localization. □ 

Appendix D 
Proof of Theorem [3] and Corollary [T] 

Proof: When the a priori knowledge of the agents' 
position is available, the derivation of EFIM, equation dZSl l 
becomes 



k£K jeAfbUMAlk} 



\nf(rkj\pk,P],Kkj) 



ln/(Kfcj|Pfc,Pj! 



ln/(P) 



Following the notations and derivations in Appendix IB -Al we 
obtain the EFIM given by (fT4l i. This completes the proof of 
Theorem [3] Note that the structure of ( fT4] i is similar to that 
of dJTT i except the additional term Hp. 

The EFIM in (fl4] i is applicable to general case. Note that 
Rfe(rfcj) in this case cannot be further simplified as that in ( |34| ) 
since we need to take expectation over the random parameter 
P in ( l32b . However, when condition ( fTSI l holds for functions 
^k] {dkj,dk]) Qfej <lly <lkj ^kj{dk],Pk), and ^kj{Kkj,Kkj), 
the expectations of those functions with respect to P can be 



iL *fej \kj + c^^k]{dkj,dkj) 



Ifej + c'Hfej {dkj , Kkj ) 1 I + r/Skj {Kkj , Kkj ) 1 ( Vkj + r/Skj {dkj , Kkj ) 



(35) 



(36) 
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replaced by the values of the functions at P. In such a case, 
the RI in (fTSl l can be written as 



where Xkj is the RII given in ( |35] | evaluated at P, and (pkj is 
the angle from pfc to pj. □ 



Appendix E 
Proof of Theorem|4] 

Proof: Consider a cooperative network with agents, 
whose overall EFIM is given by ( fT4l i. If agent A^a 
has infinite a priori position knowledge, i.e., Hp„ = 
linij2_^g^ diag {t^, P}, then we apply the notion of EFI to 
eliminate the parameter vector p^v,, in (fl4] i and have 



Je(pi,...,pjV.-l) = [Je(P)] 



2(A',-l)x2(Ar,-l) ' 



(38) 



where we have used (|39] l shown at the bottom of the page. 

Note that if we let A^^' = A^b U {A^a}, AA^' = A4\{Ara}, and 
Rfc(i-fc,JvJ = Rfc(rAr„,fe) + Rfc(rfc,jvJ for k € in dH, 
the structure of ( |38] | becomes the same as that of ( fT4l i. with 
a dimension decrease by 2. Therefore, the new RI Rj.(rA;jv, ) 
is fully utilizable, i.e., agent with infinite a priori position 
knowledge is effectively an anchor. □ 

Appendix F 
Proofs for SectionITv] 

A. Proof of Proposition |2] 

Proof: If the current coordinate system is rotated by angle 
(f) and translated by po = [ Xq ?/o then the position of 
the agent in the new coordinate system is p = p + po. 
Consequently, the EFIM for p is 



Je(p) 



9p 

T 





Je(p) 



(?P 



9p 



= UTj,(p)U^ 



Due to the cyclic property of the trace operator 
immediately find that 

T'lp) =tr{[Je(p)]-^} =tr{[Je(p)]-^} =7'(p). 



(40) 
, we 

(41) 

□ 



B. Proof of Proposition |4] 

Proof: Without loss of generality, we focus on the first 
agent. 

Lower Bound: Consider the EFIM Je(P) shown in ( l42b 
at the bottom of the page. It can be obtained from Je(P) 
by setting all C^j = for 1 < fc, j < A^^. This EFIM 
corresponds to the situation where cooperation among agents 
2 to A"a is completely ignored. One can show using elementary 
algebra that Je(P) ^ Je(P), which agrees with intuition 
since the cooperation information among agents 2 to A'a is 
not exploited. Applying the notion of EFI, we have the EFIM 
for the first agent as 

Je(Pl) = J^(Pl) 

+ [Cm-Ci,, (J^(p,) + Ci,,)-' Ci,, 

jeA4\{i} 



Since Ci j 



[cos^i.j siii^ij ] , we can express Jg(pi) as 

jL(pi) = jf(pi)+ ^L-cm' 

where ^i,, " 1 - ^i.j^l,., (Je^(P.) + Ci,,) ^ 
coefficient can be simplified as 



(43) 



q^i,,. The 



(44) 



where 



1 



1 



cos^ {'dj - 4)i,-j) H sin^ {dj ~ 



Upper Bound: Consider the EFIM J^(P) shown in ( |45] ) at 
the bottom of the next page. It can be obtained from Je(P) 
by doubling the diagonal elements C^j and setting the off- 
diagonal elements —Ckj = for 1 < A;, j < A'a. One 
can show using elementary algebra that J^(P) ^ Je(P), 
which agrees with intuition since more cooperation informa- 
tion among agents 2 to A"a is assumed in (|45] |. Applying the 



lim I ^ R7v.,(rjv„j) + ^ [RAf^(rAr_,j) + Riv.,(i"j,jvj ] 



= 



(39) 



Je(P) 



J^(Pl)+ E Ci,, -Ci,2 

jeA4\{i} 

-Ci,2 J^(p2) + Ci,2 







(42) 







J^(p^J + Ci, 
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notion of EFI and following the similar analysis leading to 
( |43] l and ( |44] |. we obtain the EFIM for agent 1 as 

J^(pi)=J^(pi)+ ^M^M' 



where 



in which 



ti.i 



1 



(46) 



^ji'Pi-j) = cos^ (dj - + sin^ (dj - 
with 11 j, fjj, and satisfying 

feeA^a\{ij} 



□ 



Appendix G 
Proof of the Scaling Laws 

Lemma 1: Let 0i's be N i.i.d. random variables with uni- 
form distribution over [0, 27r). Then, for any < e < 1, there 
exist an iVo e such that y N > Nq, 



\ N N 2 I 

Proof: First, we note that replacing with (f>i mod tt 
preserves the value of sin^(0fc ^ 4'j)- Hence, we can consider 
0i's to be i.i.d. and uniformly distributed in [0, tt). 

We order the N (f>i's, such that < < (f>(2) < ■ ■ ■ < 
0(7V) < Using order statistics ll50l . we find that the joint 
PDF of the (blip's is 



/(0(l),'/'(2),---,'/'(W)) - -]v l{O<0(i)<</.(2)<---<0(jv)<^} 



where 1 is the indicator function. From ( |48T l, the marginal 
PDF of can be derived as |501 



1 



TT^ (fc-l)!(iV-fc)! 



Now consider a large = 8A' for some integer K, and let 5 ^ 
tt/6. The function f^^j^^{x) has a maximum at a; = 7r/8, and 
is monotonically decreasing in [tt/S, tt) D [5, tt). Therefore, 
we have 



>5) <{-K-5)U,^{5) 



(49) 



Since Wxivk^oo lini^ i^) = 0' there exists ATi e N such that 
V{(I)(K)> 5] < e/4, VA: > A'l. Note also that 

P {0(7if+i) < ^ - (5} < (tt - (5) - ^) , 

and hence, for the same ATi, P{0(7k+i) < tt — (5} < e/4, 
VAT > ATi. Similar arguments show that there exists 
A'2 G N such that P {0(3x+i) < 7r/2 - (5} < e/4 and 
P {0(5K) > ^/2 + 5} < e/4, Vi^ > K2. 

Combining the above results, we have with a probability 
1 - e, 

[[0, <5], .7 = 1,..., -A, 

e <^ [Va-'J, 7r/2 + 5], j = 3A' + l,...,5-A, 

when AT > max{A'i, A'2}. Therefore, 

N N 

k=l j=k+l 

K 5K 
^E E Sill2(^(fc) -^y)) 
k=l j=ZK+l 

5K 8K 

+ E E sin2(0(fc) 

k=3K+lj=7K+l 

(K 5K 5K 8K \ 

E E 1+ E E 1 
fe=lj=3A'+l k=3K+lj=7K+l J 

= (50) 

p 

where > denotes an inequality with probability approaching 
one as K 00. Substituting A^ = 8 A', and noting that the 
summation in ( |50l ) considers only half the terms (with j > k), 
we arrive at ( |47] i. 

Moreover, the probability in ( |49] l decreases exponentially 



(48) with K, because if letting ax = fifuK^iS), 



aK+i _ (8A' + 8)(8A' + 7)---(8A: + 1) l/sV 



K^oo an (7A' + 7)(7A' + 6) • • • (7A: + l)K 6 V6, 

<1, (51) 

and hence one can see that e in (07]) decreases exponentially 
with N. □ 
Lemma 2: Let A^'s be A^ i.i.d. random variables with arbi- 
trary distribution on the support [0, Amax]- If ^ {\ < Aq} < 
e < 1/2 for some Aq e [0, Amax], then 



P{A 



(Ar/2+l) 



< 



Ao} < 



(52) 



Je"(P) 



J^(Pl)+ E Ci,,- 

— Ci_2 



— Ci_2 

J^(P2) + Ci,2+ E 2C2J 
JGA4\{1,2} 







(45) 
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where A(i) is the order statistics of such that < A^i) < 
A(2) < ■ • • < A(jv), and e ^4e (1 - e). 

Proof: Denote the probability density and distribution of 
by fx and Fx, respectively. Consider N = 2K for some 
integer K and x £ [0, A^ax] such that Fx{x) < 1/2. Using 
the order statistics, we have 

N / \ 

F\...A^)= E r^Fx{xy{l~Fx{x)f-^ 

'j=K+l ^ 
N 

< J2 2^Fx{xy{l~Fx{x)f-' 

j=K+l 



<2^{i~Fx{x)f y: 

j=K+l 

Fxix) 



Fxix) 



1-Fxix) 
< [V'^Fx{x){l-Fx{x))y , 



where the first inequality follows from (j) < 2^, the second 
inequality is due to the extension of finite summation, and the 
last inequality follows from Fx{x) < 1/2. Replacing x with 
Ao gives (|52] |. □ 



A. Proof of Theorem \5\ 

Proof: We consider first the non-cooperative case, fol- 
lowed by the cooperative case. In either case, without loss of 
generality, we focus on the first agent at position pi. 

Non-cooperative case: We will show that 7'(pi) G 
f}(l/7Vb) and 7'(pi) G 0{\/Nb), which implies that 'P(pi) G 
e(l/iVb)0 

For an ampUtude loss exponent b, signal powers decay with 
the distance following SNR(r) cx 1/r^^. We can express the 
RII from a node at distance r as 

where ro is the minimum distance between nodes determined 
by the node's physical size, rmax is the maximum distance 
between nodes determined by the fixed area associated with 
dense network setting, and random variable Z accounts for 
the large- and small-scale fading. Since Q < Z < zi for some 
zi e M+, there exists e {Q,zi) such that < zo} < 
for a given S (0, 1). Thus, the RII from the jth anchor is 
bounded as < A^in < Aij < A,nax with probability 

IPIAmin < < Amax} < 1 " , 

z,/rf. 



where Amin ^o/?'max and A 
On one hand, we have 



Je(Pl) ^ An,ax ^ ^''^ 



(53) 



'^Similar to the definition of notation 0(/(n)), the notation g(n) 6 
n(/(n)) and g(n) 6 0{f{n)) denote, respectively, that g{n) is bounded 
below by ci/(n) and above by C2/(n) with probability approaching one as 
n — > oo, for some constant ci and C2. 



By the Cauchy-Schwarz inequality, we have 

tr{[Je(pi)]-^}-tr{Je(pi)}>4. 



Since the inequality (1531 1 together with the fact that 
tr{Jr((/)i j)} = 1 imply that tr{Je(pi)} < X^axN^, we have 
that 

T'lpi) =tr{[Je(pi)]}"' >4/(A„axiVb). 

Therefore, -p(pi) e n{l/Nb). 

On the other hand, for the lower bound, we first order the 
Ab RII Ai.j's, and then the probability of A(^^/2+i) ^ -^min is 
exponentially small by Lemma |2] i.e.. 



(54) 



for some constant e G (0, 1). Let A/],' denote the set of anchors 
with RII A(j) such that j > + 1, and we have that 

A,n,„ J2 Jr('^ij) ^ Je(Pi) > > 1 - ei , (55) 
ieX' J 

where the outage probability ei decreases exponentially with 
A'b. Moreover, since 



2iVb/2 











tr < 






■j 











EfceX' EjGA/;' sin {(pi.k - <l>i,j) 

(56) 



applying Lemma [T| gives 
1 



A„ 



■ tr • 



< 



128 



> 1 - £2 , 



(57) 



for sufficiently large A^b- The inequality in ( fSSl l impHes that 

-1 ' 



Vipi) < T^tr. 

/*min 



and hence P(pi) < 128/(Aniin At,) with probability approach- 
ing one as iVb oo. Therefore, 7^(pi) G 0(1/ Nh) with 
probability 1. 

Note that since both the outage probability ei in ( fSST i and 62 
in ( fSTl i decrease exponentially with iVb, the outage probability 
e(A'b) of the scaling law in (l24l) decreases exponentially with 
iVb. 

Cooperative case: For the cooperative case, we will use the 
lower and upper approximations of the EFIM from (|22] | and 
(|23] |. The upper approximation gives 

jy(pi) = J^(pi)+ E ^".Cij 

ieA/'A{i} 
^ E -^ij Ji-(<^ij)' 

jGAAbUMAil} 

where the inequality is obtained by treating all other agents 
to be anchors, i.e., = 1 (j G A^). In this case, there are 
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equivalently A^b+^a^ 1 anchors, and similar analysis as in the 
non-cooperative case shows that P{pi) G il{l/{N\, + N^)). 

On the other hand, from the lower approximation, we have, 
with probability approaching one, that 



J^(pi)=J^(pi) 



jeA^,\{i} 



(58) 



]e^/bU^^,\{^} 

where i> > is a given lower bound on both the RII Xij (j £ 
Afb) and the effective RII (^ji^ij {j £ A/'a). From Lemma 
|2l we can find such P for the dense network setting, because 
there exist constants < ci,C2 < +oo such that Aij > ci, 
t^ij > ci, and Aj{4>i,j) < C2 with probability approaching 
one; defining £> ^ ci/(l + ci • C2) implies Ai,^ > ;> and 
> i> since ^J-^ = [ 1 + Aj((/)ij)i/ij ]^ . Applying 
Lemma [T] and |2] and following a similar line of reasoning as 
in the non-cooperative case, we find 7'(pi) G 0(l/(A'b + VVa)) 
with probability approaching one as A'b, A^a ^ +00. Thus, 
we conclude that the SPEB in cooperative networks scales as 
e(i/(A^b + iVa)). □ 

B. Proof of Theorem |6| 

Proof: Let pb denote the density of anchor nodes uni- 
formly distributed in an extended network. Consider an area 
within distance R to agent 1, then the expected number of 
anchors within that area is iVb = p\,7tR^. Following a similar 
analysis leading to ( |54| ). we can show that the effect of large- 
and small-scale fading together with path-loss on the RII can 
be bounded as ci/r^^ < X{t) < C2/r^'' for some constants 
< ci < C2 < +00, with an outage probability exponentially 
decreasing with A'b and iVa. This implies that, with probability 
approaching one, the large- and small-scale fading will not 
affect the scaling law0 and hence we can consider the RII 
from a node at distance r as 

for the analysis of the scaling laws. Since each anchor is 
uniformly distributed in the given area, the PDF of the RII 
can be written as 

1 1 

/(A) = 



b{R^-rl) 



X- 



1 



with mean 



E{A} 



6-1 



1 



1-h 



'0 

15 



and second moment 



EjA^} 



2fc-l 
1 



26-1 



{l/R^>'<X<l/r^''} I 

b=l, (59) 
< 6 < 1, 

b>l, 

b^l, (60) 
0< fe< 1. 



"it will be shown that the overall outage is dominated by the spatial 
topology for a large number of nodes, and thus we can ignore the outage 
due to fading. 



Note that N\, oc R^, we can show that the mean scales as 

'e(i/7Vb), b>i, 

E{A} e e(logiVb/A^b), b = l, (61) 
Qil/N^), 0<6<1, 

and the variance always scales as 

Var{A} e e(l/iVb) • (62) 

When > 1, it follows that, for fixed densities of anchors 
and agents, tr{Je(pi)} G ©(1) with probability approaching 
one as A'b +00, which implies that 7'(pi) G ©(I)- 

We will show that when b = I, the 'P(pi) scales as 
8(l/log A'b) and 8(l/log(A'b + A'a)) for the non-cooperative 
case and cooperative case, respectively. Using a similar ar- 
gument, we can easily show that for < 6 < 1 the SPEB 
scales as Q{l/N^-^) and e(l/(A'b + N.^)''^^) for the non- 
cooperative case and cooperative case, respectively. 

Non-cooperative case (b = 1): We introduce a random 
variable Yn, = Ejga^ ^ij/ log(A^b)- From dSB and 
we have 

lim E{Yn,}^C, 

A'b— > 00 

for some constant C, and 

lim E{|y^,-Cp} 

A'b— s-oo 

= lim Var{rjvJ + lim |E{y^J - Cf 

A'b— >co A'b— i-oo 

+ lim 2{E{Yn,}~C)-E{Yn,~E{Yn,}} 

A'b— > 00 

= 0. 

This implies that J^jeAfb ^'^■i ^^^^^^ ©(log Aj,) with prob- 
ability approaching one, and hence tr{Je(pi)} G 6 (log Ab). 
Using a similar analysis as in Appendix IG-AI we can show 
that P(pi) e 17(l/logA'b). 

For the upper bound, using the same argument as in Lemma 
[T] we can show that with probability approaching one, there 
are Ab/8 anchors with angle 4>k & [0, tt/6 ] and Ab/8 anchors 
with angle (pk G [■''■/3,7r/2] to the agent. We denote these two 
disjoint sets of anchors by JVi and A/2, and define 

Je(Pl) = E Jr((?!)l,j), 

je^flU^f2 

and 



j:(pi) " E -^1^^ + E ^1^^ 

Then, we have 
tr{[Je(pi)]-'} <tr|[je(pi)]"'| <tr|[j:(pi) 



(63) 

where the first inequality comes from JVi U Af2 ^ Afb, and the 
second inequality is due to the fact that the SPEB increases 
if we set (pij = 7r/6 for j G Afi and 0ij — tt/3 for 
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j G 7V2E3 Since both EjeM -^ij and Eje/Zs '^i-J ^'^^^^ 
0(logA'b), 'P(p) e 0(l/log7Vb) with probability approaching 
one. Therefore, the SPEB in non-cooperative extended net- 
works scales as 9(l/logA'b). 

We finally check the probability of outage, i.e., X^jeA/i -^i.i 
is not in 0(logA'b). For a fixed large A'b, the distribution 
of X]jeA4 ^i j /v'^ '^aii be approximated as the normal 
distribution 7V(log iVb///Vb, l/A'b), and henc^ 



E 



Aij -logiVb 



2g 



1 



> ^logA^b 



loaiVh 



(64) 



where Q( ) is the tail probability function of standard normal 
distribution. Approximations and bounds for the tail proba- 
biUty function can be found in ||5TI - If53l . Moreover, when 
< & < 1 a similar argument leads to 



E 

fceA^b 



At - iV. 



1-6 



1-b 



exp 



2-26 



(65) 



Cooperative case (b = 1): The cooperative case can be 
proved similar to the above non-cooperative case in conjunc- 
tion with the cooperative case of Theorem |5] It turns out that 
the SPEB can be shown to scale as 0(l/log(A'b + N-^)) when 
all other agents are considered to be anchors. We can also 
show that, with probability approaching one, the SPEB scales 
as 0(l/log(A^b + ^a)), using the lower approximation of the 
EFIM, and an argument similar to ( |63] |. □ 
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